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Abstract. In an ongoing effort to identify and study high-mass protostellar candidates we have observed in various 
tracers a sample of 235 sources selected from the IRAS Point Source Catalog, mostly with 6 < —30°, with the 
SESF antenna at millimeter wavelengths. The sample contains 142 Low sources and 93 High, which are believed 
to be in different evolutionary stages. Both sub-samples have been studied in detail by comparing their physical 
properties and morphologies. Massive dust clumps have been detected in all but 8 regions, with usually more 
than one clump per region. The dust emission shows a variety of complex morphologies, sometimes with multiple 
l^-} | clumps forming filaments or clusters. The mean clump has a linear size of ~ 0.5 pc, a mass of ~ 320 Mq for a dust 

temperature Td = 30 K, an H2 density of 9.5 x 10 5 cm -3 , and a surface density of 0.4 gem -2 . The median values 
are 0.4 pc, 102 Mq, 4 x 10 4 cm -3 , and 0.14 gem -2 , respectively. The mean value of the luminosity-to-mass ratio, 
L/M ~ 99 Lq/Mq, suggests that the sources are in a young, pre-ultracompact Hn phase. We have compared the 
millimeter continuum maps with images of the mid-IR MSX emission, and have discovered 95 massive millimeter 
clumps non-MSX emitters, either diffuse or point-like, that are potential prestellar or precluster cores. The physical 
properties of these clumps are similar to those of the others, apart from the mass that is ~ 3 times lower than 
for clumps with MSX counterpart. Such a difference could be due to the potential prestellar clumps having a 
lower dust temperature. The mass spectrum of the clumps with masses above M ~ 100 Mq is best fitted with a 
power-law dN/dM oc M~ a with a = 2.1, consistent with the Salpeter 119551 stellar IMF, with a = 2.35. On the 
other hand, the mass function of clumps with masses 10 Mq < M < 120 Mq is better fitted with a power law of 
slope a = 1.5, more consistent with the mass function of molecular clouds derived from gas observations. 
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1. Introduction 

Massive stars (M > 8 Mq) play a crucial role in the 
appearance and evolution of galaxies. They are respon- 
sible for the production of heavy elements and influence 
the interstellar medium through energetic winds and su- 
pernovae. Despite their importance, the understanding of 
massive star formation has remained significantly behind 
that of their lower-mass counterparts, for which much ob- 
servational and theoretical work has already been done. 
This situation has changed in recent years, when the for- 
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mation of high-mass stars has been gaining increasing in- 
terest, with the attention gradually shifting from the study 
of clouds associated with ultracompact UC H11 regions, to 
those where only luminous IRAS sources without radio 
continuum emission were detected. This is equivalent to 
approaching the earliest stages of high-mass star forma- 
tion, when most of the luminosity of the newly formed 
(proto)star is derived from the release of gravitational en- 
ergy. 

Palla et al. (11991|l have used the IRAS Point Source 
Catalog (IRAS-PSC) to select plausible candidates of mas- 
sive (proto)stars with 8 > —30°. The basic criteria used by 
these authors to select the targets are based on the IRAS 
colours and the lack of association with H11 regions: the 
former constraint is derived from the study of Richards 
et al. l|1987f) . who used IRAS colours to identify compact 
molecular clouds; the latter is aimed at biasing the sam- 
ple towards high-mass Young Stellar Objects (YSOs) in a 
very early phase of their evolution, when an Hn region has 
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not yet developed. The resulting sample was split into two 
sub-samples, this time using the IRAS colour selection cri- 
teria of Wood & Churchwell 1|1989[1 for identifying UC Hn 
regions: [25 - 12] > 0.57 and [60 - 12] > 1.30. The two 
sub-samples of compact molecular clouds satisfying and 
non-satisfying the Wood & Churchwell criteria have been 
called High and Low, respectively. The selected sources 
also verify that there are no upper limits for their fluxes 
at 25, 60, and 100 ^m, and that F 60Atm > 100 Jy. Palla et 
al. I|1991f) have searched for H2O maser emission associ- 
ated with the sources and found a lower association rate 
for the Low sources that has been interpreted as an indi- 
cation that the Low sources are in an earlier evolutionary 
phase than the High sources. Both sub-samples have been 
systematically observed in different continuum and molec- 
ular line tracers, from centimeter to near-infrared wave- 
lengths (Molinari et al. HW ITTm Il998hl 1217)01 mm 
Brand et al. IMfll Zhang et al. IMfll EFJ3)- The scope 
was to derive the physical properties of the two groups 
and confirm that they are in different evolutionary stages. 
In particular, the Low sub-sample would contain a fraction 
of young sources that are not yet zero-age main sequence 
(ZAMS) massive stars. The main findings of this study 
have been thoroughly discussed by Brand et al. (2001) 
and can be summarized as follow: 

— High and Low sources have luminosities typical of high- 
mass YSOs (L > 10 3 i o ); 

— the percentage of Low sources not associated with 
UC Hn regions is higher (76%; see Molinari et 
al. I1998a|l than that of High sources (57%); 

— a large fraction of Low sources have dust temperatures 
of - 30 K (Molinari et al. ESffiiUjl much lower than 
those measured towards "hot cores" (^ 100 K), and 
are hence most likely high-mass objects too young to 
have yet developed an UC Hn region. 

In order to extend our systematic search for massive 
protostellar candidates to the entire sky, we have selected 
new targets in the southern hemisphere, namely High and 
Low sources with S < —30°. This is the second of a series 
of papers aimed to conduct in the southern hemisphere 
the same kind of investigation carried out for sources with 
S > —30°. In the first paper (Fontani et al. |2005 ) we have 
discussed the results of the molecular line survey towards 
the Low sources with S < —30°. The main results of that 
study are that there is a tight association of the sources 
with dense gas, and that the physical properties of the Low 
sources, such as linewidths, and distribution of the NRAO 
VLA Sky Survey NVSS-to-IRAS flux ratio, are compara- 
ble to those found by Sridharan et al. ( 2002 ) for a sample 
of High-like sources when the luminosity of the sources is 
L < 10 5 L©. The mass-luminosity ratios are also similar 
to those found by Sridharan et al. (2002J) but lower than 
the ratio found for a sample of UC Hn regions, supporting 
the idea that our Low sources, as well as those High-like of 
Sridharan et al. (2002), are younger than UC Hil regions. 
In the present paper we discuss the main findings of the 
millimeter continuum survey carried out towards a sam- 



ple of High and Low sources with S < —30°, plus a few 
additional sources with —30° < S < 40°. In particular we 
study the morphology of the dust emission and derive the 
physical properties of the sources, and compare them with 
those reported by other surveys in the literature. In addi- 
tion, we also compare the properties of both High and Low 
sub-samples, as well as those of millimeter sources associ- 
ated with mid-infrared sources from the Midcourse Space 
Experiment (MSX 1 ) Point Source Catalog and those that 
are not. Finally, we present a study on the mass spectrum 
of the observed sources. 



2. Sample 

The first step in extending the search for massive YSOs 
towards the southern hemisphere was to select a sample 
of possible candidates from the IRAS-PS C following the 
Palla et al. I|1991|l criteria. Taking into account the interest 
of studying the earliest stages of high-mass star formation, 
we first selected a sample of Low sources with 8 < —30°, 
which was observed in C 17 and/or CS and the results 
are presented in the first paper by Fontani et al. |2005). 
Out of the 131 sources of this sample, 125 were then ob- 
served in the continuum at millimeter wavelengths (the 
results of that survey are presented in this paper). In order 
to conduct a comparative study of their properties with 
those of possibly more evolved sources, we also observed 
in the continuum a comparable sample of High sources 
with S < —30°, which had been previously detected in CS 
by Bronfman et al. I|1996(l . In addition we also observed a 
number of Low and High sources with 6 > —30°. The re- 
sulting sample observed in the millimeter continuum con- 
tains a total of 235 sources: 142 of them Low and 93 High. 
Table^shows the position, distance, and luminosity of the 
sources in the sample. As already mentioned, almost all 
(89%) the sources in our sample have 6 < —30°, whereas 
the Palla et al. l(B91) sample contained only objects with 
5 > —30°. Therefore, one may reasonably expect that our 
sample contains a higher contamination by Hn or UC Hn 
regions than that of Palla et al. 1)1991(1 . because radio con- 
tinuum surveys of Hn regions with S < —30° are less nu- 
merous than and not as complete as those with S > —30°. 

All the High sources in the sample have been previously 
detected in CS by Bronfman et al. (1996), and those with 

5 > -30°, with the exception of IRAS 18198-1429, also 
in NH3 by Molinari et al. H1996J) . The Low sources with 

6 < -30°, with the exception of IRAS 15579-5347, have 
been observed in CS by Fontani et al. (2003, an d some 
of them have also been observed in C 17 0, and those with 
S > -30° in NH 3 by Molinari et al. (fT5M| . Fifteen of the 
Low sources with S < —30° were not detected either in 
C 17 or CS (Fontani et al. l2UU5|l . 



MSX images have been taken from the on-line MSX 
database http: //www . ipac . caltech.edu/ ipac/msx/msx .html 
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3. Observations 

The 1.2-mm continuum observations were carried out 
with the 37-channel bolometer array SIMBA (SEST 
Imaging Bolometer Array) at the SEST (Swedish-ESO 
Submillimetre Telescope), on July 16-20, 2002 and July 
9-13, 2003. 

Maps of 900" x 400" (azimuth x elevation) around 
each IRAS source in Table Q] were obtained, with a scan 
rate of 80"/s, and a separation of 8" in elevation between 
scans. Bigger maps were also obtained by using a mosaic- 
ing technique towards sources 16153—5016, 17225—3426, 
and 18014—2428 (here and in the following we will refer 
to the sources without mentioning "IRAS" in the name). 
For the region surrounding source 16428—4109, the cen- 
ter of the observations was shifted ~ 4.'7 northeast from 
the nominal IRAS point source position: after taking a 
first map centered at the IRAS position, no millimeter 
emission was detected at the source nominal position but 
there clearly was emission at the edge of the map, so 
we decided to shift the center towards the northeast. We 
have checked the High RESolution (HIRES) IRAS im- 
ages and have found that there is an infrared source at 
the nominal IRAS point source position, which indicates 
that the position given in the IRAS-PSC is correct. In our 
study we have considered as associated with the IRAS 
source all the millimeter clumps located at < 90" from 
the nominal IRAS point source position. Therefore, the 
clumps detected towards 16428—4109 are not associated 
with the IRAS source. The total integration time per map 
was about 15 minutes, and the typical noise level in the 
maps is 25-40 mJy beam -1 . Atmospheric opacity was de- 
termined from skydips, which were taken every 2 hours; 
values at zenith ranged between 0.21 and 0.50 (during the 
2002 observations) and 0.13 and 0.30 (in 2003). The data 
were calibrated using observations of Uranus, made once 
or twice per day; the conversion factor ranged between 58 
and 75 mJy/count in 2002, and between 50 and 69 in 2003. 
The calibration uncertainty is about 15%. The pointing of 
the SEST was determined to be accurate within a few arc- 
sec by observing the strong continuum sources rj Carinae 
or Centaurus A every 2 hours. The HPBW is ~ 24". 

All data were reduced with the program MOPSI, writ- 
ten by R. Zylka (IRAM, Grenoble), according to the 
instructions given in the SIMBA Observers Handbook 
(2003). See Chini et al. (2003) for a detailed description 
of the data reduction method. 

4. Results 

4.1. Kinematic distances 

Table ^ lists the kinematic distances, d, to the IRAS 
sources. The distances have been estimated, using the ro- 
tation curve of Brand & Blitz l|1993[) . from the CS line 
velocity given by Bronfman et al. (1996) for all the High 
sources, from the CS line velocity given by Fontani et 
al. (jSUOljjl for those Low with 5 < -30°, and from the NH 3 
line velocity given by Molinari et al. (1996) for the Low 



sources with 8 > —30°. This method is valid for galacto- 
centric distances between 2 and 25 kpc. For sources inside 
the solar circle, there are two solutions for the kinematic 
distance, near and far. In some cases this ambiguity can 
be solved, for example, when the height of the source from 
the galactic plane is more than 200 pc (i.e. roughly twice 
the scale height of the molecular disk) at the far distance, 
or when the near distance is too small for a massive star 
forming region. This is the case for sources 14394—6004 
and 17040—3959, for which d ncar is < 100 pc, and there- 
fore the far distance was adopted. There are many sources 
for which it was not possible to solve the distance ambigu- 
ity: for these in Table^we report both distance estimates. 
Note that in the following to derive the physical parameters 
of the sources and in case of unsolved distance ambiguity, 
the near distance is adopted. The distances to the IRAS 
sources are between 130 pc and 27.1 kpc, with an average 
value of ~ 4.5 kpc, and a median value of ~ 3.8 kpc. 

We have not been able to derive an estimate of the dis- 
tance for sources 08488-4457, 10088-5730, 10156-5804, 
10575-5844, 11431-6516, 11476-6435, 12434-6355, 
13078-6247, 13558-6159, 14198-6115, 14412-5948, 
15506-5325, 15579-5347, 16204-4943, 16581-4212, 
17140-3747, 17230-3531, 17231-3520, 17242-3513, 
17352-3153, 17410-3019, 17425-3017. There are three 
possible reasons for that: no estimate of the systemic ve- 
locity of the region was available; the corresponding dis- 
tances estimates are outside the galactocentric interval 2- 
25 kpc; the height from the galactic plane is too large 
(> 200 pc) for near and far distances. 

Most of the regions observed contain more than one 
millimeter clump (see Section f4.3|l . and it is likely that 
all of them belong to the same star forming region as 
the IRAS source. Therefore, when deriving their physi- 
cal parameters, we have made the assumption that all the 
clumps are located at the same kinematic distance as the 
one estimated for the IRAS source. 

4.2. Luminosities 

The bolometric luminosities in Table were calculated 
by integrating the IRAS flux densities. The contribution 
from longer wavelengths was taken into account by ex- 
trapolating according to a black-body function that peaks 
at 100 /im and has the same flux density as the source at 
this wavelength. The distribution of luminosities is shown 
in Fig. H The average luminosity is ~ 6.7 x 10 4 L Q and 
the median is ~ 1.6 x 10 4 L Q . These luminosities are ~ 3.5 
times lower than the average and median values derived by 
Faundcz et al. (2004) for their sample of southern sources. 
This could be due to the fact that their sample has been 
selected from the survey of Bronfman et al. (1996), which 
contains high-mass very YSOs but it is also contaminated 
by more evolved sources such as Hn regions (or equiv- 
alently more massive stars), which are expected to be 
brighter at far-infrared wavelengths. 
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Fig. 1. Histogram of the luminosity distribution calcu- 
lated from the IRAS flux densities of the sources. 



4.3. Identification of the clumps 

Massive dust chimps have been detected in all 
but 7 regions, 10102-5706, 11476-6435, 14198-6115, 
15571-5218, 16403-4614, 16417-4445, and 18024-2231, 
of our sample. As already mentioned in Section for 
16428—4109, the clumps are displaced by ~ 4'7 from the 
position of the IRAS source. The millimeter maps of our 
survey are shown in Fig- El Figure |3| shows the two largest 
areas mapped in our survey, that is, the mosaiced maps 
towards sources 17225-3426 and 18014-2428. As can be 
seen in the maps, there is usually more than one clump per 
region. In some cases the large number of clumps and the 
extended emission detected in the region make it very dif- 
ficult to separate them from each other by eye. Therefore, 
in order to identify the millimeter clumps and their prop- 
erties adopting a more objective criterion, we have used a 
two-dimensional variation of the clump-finding algorithm 
Clumpfind developed by Williams et al. |1994f . The three- 
dimensional version of the algorithm was originally de- 
signed to be applied to spectral line datacubes, and the 
two-dimensional version is a simple modification of it. The 
algorithm works by effectively contouring the data at a 
multiple of the rms noise of the map, then searching for 
peaks of emission to locate the clumps, and finally fol- 
lowing the clump profile down to lower intensities. The 
contouring levels have to be chosen by hand, which means 
that the clump-finding procedure is not completely au- 
tomated, and therefore one can introduce biases into the 
results. Clumpfind does not a priori require any particu- 
lar shape of the clump profile, as some other clumps find- 
ing algorithms do, and one of its disadvantages is that it 
misses low-mass clumps that lie below the lowest contour. 
This could flatten the low-mass end of the mass spectrum 
of the regions. In our case, we set the lowest contour level 
to 3cr, and then increased the contouring by steps of 3a. 

The clump-finding procedure calculates the peak po- 
sition, the full width at half maximum (FWHM) not cor- 
rected for beam size for the x-axis, FWHM X , and for the y- 
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Fig. 3. Overlay of the 1.2 mm continuum emission (con- 
tours) obtained with SIMBA at the SEST antenna, 
on the MSX emission at 21 /im (Band E) (image) 
towards the sources IRAS 17225—3426 (upper panel) 
and IRAS 18014-2428 (bottom panel). The contour lev- 
els range from 0.11 (3<r) to 1.67 Jy beam -1 in steps 
of 0.22 Jy beam -1 (upper panel) and from 0.09 (3a) 
to 3.20 Jy beam -1 in steps of 0.17 Jy beam -1 (bot- 
tom panel). The star marks the position of the IRAS 
sources. The coordinates are in J2000 epoch. In the bot- 
tom panel, the two IRAS sources correspond to Mol 36 
(IRAS 18014-2428) and Mol 37 (IRAS 18018-2426) from 
Molinari et al. (1996). The dashed polygon shows the lim- 
its of the area mapped with SIMBA. 



axis, FWHM y , and the total flux density integrated within 
the clump boundary, that is within the lowest contour 
level. Table gives the offset positions in arcsec with 
respect to the nominal IRAS point source position for 
each clump, the angular diameter, which has been calcu- 
lated as the deconvolved geometric mean of FWHM X and 
FWHM y , the linear diameter, the total flux density, the 
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4.4. Morphologies of the clumps 
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Fig. 4. Histograms of some parameters of the clumps de- 
tected in the continuum emission at 1.2 mm: (a) kinematic 
distance, (b) number of clumps per map, (c) radius of the 
clumps, (d) mass of the clumps, (e) H2 volume density, 
and (f) H2 surface density. Note that in the regions for 
which mosaiced observations were carried out, we have 
included all the clumps in the field when computing the 
histograms, and not only the sources in a field 900" x 400", 
like for the other regions. 



mass, and the density. The last column indicates whether 
the clump is associated with MSX emission, either point- 
like or diffuse, or not. Figures 0Ji an d Eb show the his- 
tograms with the kinematic distance of the clumps and 
the number of clumps per region. The mean and me- 
dian values of the distance of the clumps are 3.9 kpc and 
3.4 kpc, respectively, while those of the number of sources 
per region are 2.8 and 2.0, respectively. Note that the val- 
ues of the mean and median distance are slightly lower 
when taking into account all the clumps in the regions in- 
stead of taking only into account the central IRAS source 
(see Section 14.1(1 . The maximum number of clumps have 
been detected in the large mosaiced maps around sources 
17225-3426 and 18014-2428, for which we have found 27 
and 22 clumps, respectively (see Fig- - 



The 1.2 mm maps in Fig. [21 show that the dust emission 
towards massive star forming regions presents a variety of 
complex morphologies. The emission, which usually peaks 
at or near the IRAS position, is very often associated with 
multiple compact embedded objects and extended emis- 
sion. Sometimes the emission is very centrally peaked, 
with a single massive clump associated with the IRAS 
source, e.g. 12063-6259, 15454-5335, and 18144-1723, 
sometimes with more than one massive clump clus- 
tered towards the central position, e.g. 12127—6244, 
13333—6234, and 17149—3916, or with smaller and fainter 
clumps clustered or located nearby, e.g. 15015—5720, 
15519-5430, 17082-4114, and 17118-3090. Sometimes 
the emission shows multiple clumps located throughout 
the field, which can be separated from each other, e.g. 
16085-5138, 16093-5128, 16153-5016, and 17242-3513, 
or linked in a chain of clumps and elongated struc- 
tures, e.g. 10184-5748, 13039-6108, 14000-6104, and 
16164—4929, suggesting star formation in a sheet or a fil- 
ament (e.g. Larson fl985j) . Particularly nice examples of 
this latter phenomenon are seen in the mosaiced maps of 
17225-3426 and 18014-2428 (Fig. EJ) that show strings 
of emission clumps and elongated and clustered struc- 
tures. Sometimes the emission towards the IRAS posi- 
tion is very faint (e.g. 08488-4457, 10277-5730, and 
15464—5445), or no millimeter emission is detected at 
all. In particular, the number of IRAS sources not de- 
tected at 1.2 mm is 12, taking into account 10545—6244, 
16153-5016, 16428-4109, and 17156-3607 as well, for 
which although there are millimeter clumps detected in 
the region, none of them is associated (i.e. within 90") 
with the IRAS source. All of them are Low sources. As 
can be seen in Table 3 of Fontani et al. (2005J), 5 of them, 
16403-4614, 16417-4445, 16153-5016, 16428-4109, and 
17156-3607 have been detected in CS and C 17 0, 3 of 
them, 10102-5706, 15571-5218, and 10545-6244 have 
been detected in CS but not in C 17 0, and 3 of them, 
11476-6435, 14198-6115, 14412-5948, have not been de- 
tected in CS but not been observed in C 17 0. Finally, one 
of the sources, 18024—2231, has been detected in NH 3 
by Molinari et al. ©96). On the other hand, there are 
6 Low sources, 08488-4457, 10088-5730, 12434-6355, 
15506-5325, 16204-4943, and 16581-4212, that have 
been detected in the millimeter continuum, although the 
emission is faint, but not in CS (Fontani et al. 2005), and 
one source, 15579—5347, that has not been observed in CS 
and has been clearly detected in the continuum. 

4.5. Linear Diameters 

The deconvolved linear diameters of the clumps have been 
computed from their angular diameters. As mentioned in 
Sect. 14.31 both linear and angular diameters are listed in 
Table El The mean value is 0.5 pc for those sources that 
have been resolved, which is in agreement with the average 
value of 0.6 pc found by Williams et al. (2004) for the 
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sources of the Sridharan/Beuther sample (Sridharan et 
al.|2D02; Beuther et al. I2002[) . and quite smaller than the 
value of 0.8 pc found by Faundez et al. l)2004[l . or 1 pc 
found by Hill et al. (2005) for their samples of southern 
sources, although it is not clear whether the sizes given by 
these authors have been deconvolved or not. The median 
value for the clumps in our sample is 0.4 pc. Figure 
shows the histogram of the radius of the clumps. 

We have searched for possible asymmetries in the 
clumps by calculating FWHM x /FWHM y . The mean and 
median values obtained, 1.04 and 0.96, respectively, indi- 
cate that the clumps are quite symmetric. 

4.6. Masses and densities 

The masses of the clumps given in Table |21 have been es- 
timated assuming that the dust emission is optically thin, 
by using 

g S v d 2 

■'"Clump — r> frri \ 1 X 1 ) 

where S v is the flux density, d is the distance to the source, 
n v is the dust mass opacity coefficient, g is the gas-to-dust 
ratio, and B U {T^) is the Planck function for a blackbody 
of dust temperature Td, all measured at v = 250 GHz. We 
adopted K250 = 1 cm 2 g _1 (Ossenkopf & Henning^94), 
g = 100, and Td = 30 K for all the sources. Estimates of 
Td have been obtained by fitting grey-bodies to the spec- 
tral energy distribution (SED) of those IRAS sources that 
have only one millimeter clump associated with them. As 
already done by Fontani et al. ( 20Q2|) f° r the Low sources 
in our sample, we have neglected the IRAS fluxes at 12 
and 25 /j,m in the fit. The reason for this choice is that 
two components are known to be present in the SEDs of 
luminous YSOs (see Sridharan et al. 2002 and Beuther et 
al. l2002f) : one associated with compact, hot gas and dom- 
inating the 12 and 25 fim fluxes; the other due to parsec- 
scale, colder material, contributing to the 60 and 100 /im 
emission. The latter component is the one of interest to 
us, because we want to estimate the temperature of the 
parcsec-scale clumps mapped at 1.2 mm. Adopting a dust 
absorption coefficient proportional to v 2 , we have found 
best fits with mean temperatures of 28 K for both Low 
and High sources. Therefore, taking into account that for 
most of our clumps there are not enough measurements at 
different wavelengths to properly fit their SEDs (for some 
of them the 1.2 mm is the only one available), we have 
decided to adopt Td = 30 K for all the clumps. A simi- 
lar value, Td = 32 ± 5 K, has been found by Molinari et 
al. ( 2000) for a sample of 30 luminous Low sources in the 
northern hemisphere. Furthermore, Faundez et al. (2004), 
who have fitted the SED of a sample of southern hemi- 
sphere sources similar to the High sources in our sample 
with grey-bodies with two components, have found an av- 
erage value Td ~ 32 K for the colder component. 

Figure 0Ji shows the histogram of the distribution of 
masses of the clumps. The mean value of the clump mass 
is 320 Mq, although the median mass is much lower, 




12 3 
Log [L(L )/M(M G )] 




12 3 

Log [L(L )/M(M )] 

Fig. 5. (a) Histograms of the distribution of the 
luminosity-to-mass (L/M) ratio for all the sources, where 
M is the sum of the masses of clumps located < 90" 
from the nominal IRAS point source position, (b) Same 
as above, for the Low (solid line) and the High sources 
(dashed line). 



102 Mq. These values are in agreement with those of 
330 Mq and 143 Mq, respectively, found by Williams 
et al. I|2004[l for the 68 high-mass protostellar candidates 



of the Sridharan/Beuther sample (Sridharan et al. IzT 
Beuther et al. 12002(1 assuming the near kinematic dis- 
tance. Williams et al. (2004) have derived their clump 
masses using the dust opacity coefficients of Ossenkopf 
& Henning (fT^Mjl . and T d ranging from 30 to 60 K. Note 
that the average value of 5.0 x 10 3 Mq given by Faundez et 
al. I)2004[l refers to the total mass in each region; i.e., the 
sum of the masses of all the clumps in that region. If we 
take into account all the clumps in each region the average 
mass is 955 Mq, which is still five times lower than the 
one derived by Faundez et al. ( 2004 ). This could be due to 
the fact that their sample contains more massive stars, as 
already suggested by the bolometric luminosities of their 
targeted sources (see Section FOI) . Hill et al. f2Q05| in their 
recent survey of massive star-forming regions harbouring 
methanol masers and/or radio continuum sources have re- 
ported an average mass for their sample of 1.5 x 10 3 Mq, 
and a median value of 1.0 x 10 3 Mq, for a dust tempera- 
ture of 20 K. The average mass would be ~ 0.9 x 10 3 M Q 
and the median ~ 0.6 x 10 3 Mq, for a dust temperature 
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Log [L/LJ 

Fig. 6. Histogram of the luminosity distribution calcu- 
lated from the IRAS flux densities, for the Low (solid line) 
and the High sources (dashed line). 

of 30 K, which is the one that we have used for our esti- 
mates. Such values are still much higher than the ones of 
our sample. This, again, could indicate that the YSOs in 
their sample arc intrinsically more massive. Unfortunately 
these authors do not report the bolometric luminosity of 
the sources, and we cannot corroborate this hypothesis. 

Figures 0^ and 0f show the histograms of the H2 vol- 
ume density and surface density of the clumps. These have 
been derived assuming that the clumps have spherical 
symmetry and a mean molecular weight of /1 = 2.29toh- 
The average values are 9.5 x 10 5 cm -3 and 0.4 gem -2 , 
respectively, and the median values are 4 x 10 cm~ 3 and 
0.14 gem -2 , respectively. 

5. Discussion 

5.1. The Luminosity-Mass relation 

An important parameter for establishing the age of a 
core is the luminosity-to-mass ratio, L/M. This ratio is 
expected to increase with time as more and more gas 
is converted into stars during the star formation pro- 
cess and the embedded cluster becomes more luminous. 
Figures [5^ and [SJd show respectively the L /M ratio for 
all the sources in our sample, and for the Low and the 
High sub-samples separately. The luminosity has been de- 
rived from the IRAS flux densities (see Sect. I4.2J) . and 
the mass is the sum of the masses of the clumps located 
< 90" from the nominal IRAS point source position (see 
Table as we have considered only these clumps as asso- 
ciated with the IRAS source. The mean value of the dis- 
tribution is 99 Lq/Mq for all the sources, and 98 Lq/Mq 
and 101 Lq/Mq for the Low and the High sub-samples, 
respectively. These values are in agreement with the value 



of 100 Lq/Mq found by Sridharan et al. pHl for the 
sources of the Sridharan/Beuther sample, when rescaled 
to the Ossenkopf & Henning (1994) opacity of 1 cm 2 g _1 
used by us to derive the masses. Sridharan et al. 1|2002[) re- 
port an average value of 0.05 M / L Q , based on the masses 
derived by Beuther et al. l2002J) . which were estimated us- 
ing the opacity of Hildebrand ({198311 . This opacity is ~ 2.5 
times smaller than the Ossenkopf & Henning opacity at 
250 GHz. However, note that the masses estimated by 
Beuther et al. B0Q2), as recently reported by these au- 
thors, are wrong by a factor of 2 (see Beuther et al. 2005 ) ; 
that is, they should be a factor of 2 lower. Therefore, the 
average value of the M/L ratio for their sample is actually 
0.025 Mq/Lq, or 40 Lq/Mq, which rescaled to the opac- 
ity of Ossenkopf & Henning is 100 Lq/Mq. These values 
are slightly higher than the average value of 71 Lq/Mq 
reported by Faiindez et al. |2004). However, as already 
mentioned, these authors have used the sum of the masses 
of all the clumps in each region in their calculations. If we 
sum the mass of all the clumps in the regions, the average 
value that we obtain is 77 Lq /Mq , a value similar to the 
one reported by those authors. 

No significant difference is seen between the mean 
value of L/M for the Low and the High sources. The same 
conclusion has been found by Fontani et al. ( 2005) when 
comparing the Low sources of our sample with the High 
sources of the sample of Sridharan/Beuther (Sridharan et 
al. EEEI Beuther et al. I2002j> . These latter authors have 
compared the L/M values of the sources in their sam- 
ple with those of known UC H11 regions with the same 
masses (Hunter IT§9T1 Hunter at al. E^UUUjl and have found 
a L/M ratio 5 times higher for the latter. Note that all 
these masses have been wrongly calculated by these au- 
thors with a factor 2 higher (see Beuther et al. |5005). 
However, this error do not change the conclusions derived 
from their comparison, as both are affected by the same 
factor. Mueller et al. <|2002(l and Faiindez et al. 1)2004[> 
have found a similar result, although the difference be- 
tween cores associated with UC Hn regions and those not 
associated reported by these authors is not statistically 
significant. Sridharan et al. (2002) have suggested that 
the difference in L/M ratio occurs because the sources in 
their sample are in a younger pre-UC Hn phase, and that 
L/M increases as the cores evolve and develop UC Hn re- 
gions. Taking into account that the sources in our sample 
have L/M ratios similar to those of the Sridharan/Beuther 
sample (see Sect. 14.611 . we conclude that they too are 
younger than UC Hn regions. 



5.2. Low versus High 

One of the goals of this work is to carry out a comparative 
study of two sub-samples of massive YSOs, Low and High 
believed to be in different evolutionary stages. To do this 
one should check whether the two groups have different 
physical properties that could confirm a different evolu- 
tionary phase. Figure shows the histogram of the lumi- 
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Fig. 7. Histograms of (a) kinematic distance, (b) number of clumps per map, (c) radius of the clumps, (d) mass of 
the clumps, (e) H2 volume density, and (f) H2 surface density, for the Low (solid line) and the High sources (dashed 
line). Note that in the regions for which mosaiced observations were carried out, we have included all the clumps in 
the field when computing the histograms, and not only the sources in a field 900" x 400", like for the other regions. 



nosity distribution for both sub-samples, and Fig.[7|shows 
comparative histograms of some other physical properties, 
such as the number of clumps in each region, the size of the 
clumps, their mass, H2 volume density, and surface den- 
sity. As can be seen in these figures, both sub-samples Low 
and High are quite similar and do not show significant dif- 
ferences. This is the same conclusion reached by Fontani 
et al. (2005) when comparing the properties (linewidth, 
luminosity, temperature, NVSS-to-IRAS flux ratio) of the 
Low sources in our sample with the High sources of the 
Sridharan/Bcuthcr sample with luminosity L < 10 £©. 

Analyzing the properties in more detail, one can see 
that the luminosity distribution of the High sub-sample 
(Fig. EJ) shows more sources at higher luminosities; the 
mean value for the High sub-sample, 1.3 x 10 5 L©, is a 
factor of 5 higher than that for the Low sub-sample, 
2.5 x 10 4 L©. The median value, which is 2.6 x 10 4 L© for 
the High and 1.3 x 10 4 L© for the Low, is only a factor 



of 2 different. The mean number of clumps per region is 
similar for both sub-samples, 2.9 for the Low and 2.7 for 
the High, while the median number is 2.0 for both sub- 
samples. The maximum number of clumps per region are 
27 and 22 and correspond to the mosaics around the Low 
sources 17225-3426 and 18014-2428, respectively. If one 
does not take into account these two regions, the maxi- 
mum number of clumps per region are comparable, with 
12 clumps for the Low source 17141—3606, and 10 for the 
High source 16085—5138. The morphologies of the clumps 
are quite similar in both sub-samples (see Fig.(2J. Not sur- 
prisingly, the sources with fainter millimeter emission or 
no emission at all towards the IRAS position, which are 
mostly Low sources, have no C 17 or CS detected (Fontani 
et al.EDD3- 

Some differences are evident in the linear size and the 
mass of the clumps, which have mean and median values 
of 0.46 pc and 0.38 pc, and 163 M© and 64 M© for the Low 
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clumps, and 0.67 pc and 0.56 pc, and 561 M Q and 164 M Q 
for the High ones. Brand et al. l|2001fl have found an oppo- 
site result for the linear diameters of sources with S > —30, 
as they have found that clumps around Low sources are 
~ 3 times larger than those around High clumps. However, 
those linear diameters have been obtained from molecu- 
lar line observations instead of from the dust continuum 
emission. 

As explained in Sect. |4~Tj1 we have derived a tempera- 
ture of ~30 K for all the clumps, with no significant differ- 
ence between Low and High sources. We have estimated 
the masses of the clumps adopting this value. However, 
our temperature estimates are based on a very limited 
number of measurements, so that we cannot rule out the 
possibility that High sources are hotter than Low ones. As 
a matter of fact, while Faundez et al. (2004) have found 
a temperature (32 K) similar to ours for their sample of 
High-like sources, Sridharan et al. H2()()2|l determined a 
higher value (50 K) for a similar sample. If one adopts 
50 K also for our High sources, the mass estimates ob- 
tained by us should be multiplied by 0.55 and the mean 
value of the mass of the clumps would be 309 M Q , which 
is still a factor ~ 2 higher than the mean value of the Low 
sources. Therefore, the difference in mass of the clumps 
seems to be real. Regarding the other physical proper- 
ties, the H 2 volume density, and the surface density of the 
clumps, the average and median values are 1.4 x 10 6 cm -3 
and 4.1 x 10 4 cm -3 , and 0.45 gem -2 and 0.13 gem -2 for 
the Low sources, and 2.6 x 10 5 cm -3 and 3.9 x 10 4 cm -3 , 
and 0.29 gem -2 and 0.16 gem -2 for the High ones. 

We have searched for possible asymmetries in the Low 
and High clumps by calculating FWHM x /FWHM y , and 
we have found similar values for both sub-samples. The 
mean and median values are 1.04 and 0.95, respectively, 
for the Low clumps, and 1.03 and 0.97 for the High ones. 
Such values indicate that the clumps of both sub-samples 
are quite symmetric. 

The two sub-samples seem to be in the same evolution- 
ary stage, as suggested by the L/M ratio (see Sect. I5.ljl . 
Therefore, in conclusion, there is no significant difference 
in the physical parameters derived from the present mil- 
limeter continuum observations of the Low and High sub- 
samples with 6 < —30°. Both seem to be hosting high- 
mass (proto)stars with similar physical properties (M, 
nn 2 , radius, L/M) and evolutionary stage. However, in 
order to reasonably compare the results of our southern 
study of Low and High sources with those obtained in 
previous northern hemisphere (5 > —30°) studies (Palla 
et al. ITWn Molinari et al. ITTM 1217)711 Brand 

et al. I2001J) which found differences between both sub- 
samples, more observations at different wavelengths of 
line and continuum emission are still needed in the south- 
ern hemisphere. In particular, it would be very important 
to properly estimate the temperature of the sources, ei- 
ther through some molecular line observations, or though 
more infrared and (sub)millimeter observations to better 
constrain their SEDs. Such temperature estimates could 
tell us whether the Low sources are colder, and therefore 
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Fig. 8. Histogram of the 1.2 mm to 21 /an integrated flux 
density ratio for clumps associated with point-like MSX 
emission (solid line). The dashed line indicates the me- 
dian value of S mm /(1.5 Jy), where S mm is the 1.2 mm 
integrated flux density and 1.5 Jy is the MSX detection 
limit at 21 /xm, for those clumps that are not associated 
with MSX emission. Note that due to the fact that these 
clumps have not been detected at 21 /jm, the median value 
of the ratio S mm /(1.5 Jy) is a lower limit. 



younger, than the High sources, as observed in the north- 
ern hemisphere. 

5.3. MSX versus non MSX 

As seen in the previous sections, the IRAS sources in our 
sample are usually associated with massive clumps likely 
hosting high-mass (proto)stars. In an attempt to approach 
even closer to the earliest stages of high-mass star forma- 
tion, it would be of great interest to identify a core prior to 
the star formation process; namely a prestellar or, better 
called, precluster core. Such a precluster core is expected 
to have density and size similar to those with embedded 
high-mass protostars, but lower luminosity and tempera- 
ture. The bulk of its luminosity is expected to be emitted 
at millimeter and submillimeter wavelengths, with faint 
or no mid-IR or far-IR emission. One may reasonably ex- 
pect that such cores capable of forming massive stars but 
in a stage prior to the onset of star formation are located 
close to massive luminous cores containing already formed 
high-mass (proto)stars. The millimeter continuum maps 
have shown that there is usually more than one millime- 
ter clump in the region surrounding the IRAS source, and 
that not all of them are associated with it; that is, they are 
located at > 90" from the nominal IRAS point source po- 
sition. Therefore, we have searched for precluster cores in 
the surroundings of the candidate massive (proto)stars by 
cross-correlating our sample with the MSX Point Source 
Catalog, since the diameters of the clumps are compara- 
ble to the spatial resolution of 18'.'3 of the mid-IR MSX 
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Fig. 9. Histograms of (a) radius of the clumps, (b) mass 
of the clumps, (c) H2 volume density, and (d) H2 surface 
density, for clumps associated with MSX emission, either 
point-like or diffuse (solid line), and those that are not 
(dashed line). Note that in the regions for which mosaic 
observations were carried out, we have included all the 
clumps in the field when computing the histograms, and 
not only the sources in a field 900" x 400", like for the 
other regions. 



observations. The cross-correlation has been done search- 
ing for mid-IR point sources in a radius of 40" around 
each millimeter clump, in any of the four bands (8 /urn, 
12 fim, 14 /im, or 21 /im) in the MSX-PSC version 6.0. 
In case that more than one millimeter clump was found 
at < 40" from the same MSX point source, we have arbi- 
trarily chosen the nearest clump as the associated one. We 
have also compared by eye the millimeter continuum maps 
with images of the MSX emission (see Fig. |3J), which has 
allowed us to identify millimeter sources associated with 
mid-IR extended or diffuse emission. In the last column of 
Table |3 it is indicated whether the clumps are associated 
with point-like or diffuse MSX emission, or have no MSX 
emission at all. The positions in the MSX-PSC are accu- 
rate to an rms of 4"-5", and the calibration uncertainties 
are well within 13% (Egan et al. [TM5jl . 

As a result of this study we have discovered 95 mas- 
sive millimeter clumps not associated with mid-IR emis- 
sion, meaning no mid-IR point source at any of the four 
bands, nor diffuse emission, that are potential preclus- 
ter clumps (see Table 0). However, this number could 
be a lower limit because the excess of mid-IR emission 
in the observed regions produces, in some cases, confu- 
sion. Therefore, the association of some millimeter clumps 
with diffuse MSX emission may not be real. In addition, 
there are 35 clumps in the MSX-PSC detected only at 



8 /fin, whose emission could be due to contamination by 
Polycyclic Aromatic Hydrocarbons (PAHs) and not due to 
an embedded source. Until now only a very limited num- 
ber of such candidates have been found in the northern 
(Molinari et al. I1998bl Beuther et al. I3UU21 Sridharan et 
al. 2002), and southern hemisphere (Garay et al. [2004 
Hill et al. 12005(1 . Note that in some cases, such as for ex- 
ample source 17221—3619, there are MSX sources located 
at the edge of the millimeter cores, which might be either 
reflection nebulae or sources belonging to the same star- 
forming region but in an older evolutionary stage: in fact, 
more evolved, and therefore hotter sources, are stronger 
emitters at mid-IR wavelengths. Seventy out of the 95 
clumps without MSX emission are clumps in maps around 
Low sources, which is the 16% of the Low clumps, while 20 
(10%) are clumps in maps around High sources. Therefore, 
the Low sub-sample contains in percentage more potential 
precluster clumps than the High sub-sample. 

Figure |H1 shows the histogram of the 1.2 mm to 21 /im 
integrated flux density ratio for the clumps associated with 
point-like MSX emission, where the flux density at 21 /zm 
has been taken from the MSX-PSC. The dashed line in- 
dicates the median value of S mm /(1.5 Jy), where S mm is 
the 1.2 mm integrated flux density and 1.5 Jy is the typi- 
cal MSX detection limit at 21 /zm, for those clumps that 
are not associated with MSX emission. Note that because 
these clumps have not been detected at 21 /im, the detec- 
tion limit flux of 1.5 Jy is an upper limit for their emis- 
sion at 21 /im, and therefore, the median value of the ratio 
S mm /(1.5 Jy) is a lower limit. As can be seen in Fig. 
there are relatively strong millimeter clumps that have no 
MSX counterpart, which indicates that they are not as- 
sociated with mid-IR emission. Therefore, some of the 95 
millimeter clumps without MSX counterpart could be in 
fact precluster clumps, where luminous stars have not yet 
formed. 

Figure shows the comparative histograms of some 
physical properties, such as the radius, mass, H2 volume 
density, and surface density of the clumps associated with 
MSX emission, either point-like or diffuse, and those not. 
As can be seen in the distributions, the physical properties 
of the two sub-samples are similar. Due to the fact that 
the scale of the distributions is logarithmic, the properties 
would have to be very different in order to be clearly vis- 
ible in the histograms. The main difference is the mass of 
the clumps, which for non-MSX emitters is clearly lower, 
with a mean value of 96 M Q and a median of 33 M Q 
compared with the mean and median values of 336 Mq 
and 106 Mq, respectively, of the clumps with MSX coun- 
terparts. However, the dust temperatures of the clumps 
without mid-IR emission could be significantly lower than 
that of the clumps with embedded massive (proto)stars, 
which is Td ~ 30 K. In fact, Garay et al. (200]]) have found 
Td < 17 K for those clumps not associated with MSX 
emission in their sample. Assuming a dust temperature of 
15 K, the masses of the clumps should be corrected by a 
factor ~ 2.5, and therefore, the masses of both MSX and 
non-MSX emitters would be more comparable. Regarding 
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Fig. 10. Left top panel: The mass spectrum of the 1.2 mm clumps detected at a distance < 6 kpc. The solid line 
represents the Salpeter IMF, dN/dM oc M~ 2 - 35 ; the dotted line is a —2.1 power law, obtained from the least square 
fit to the data, and the dashed line is a —1.7 power law. The vertical dot-dashed line indicates the completeness limit 
at 6 kpc. Right top panel: The normalized cumulative mass distribution of clumps with masses above the completeness 
limit at 6 kpc. The solid, and dashed lines are the same as in the left panel, and the dotted line is a —1.9 power law, 
obtained from the least square fit to the data. Left bottom panel: Same as above for clumps detected at a distance 
< 2 kpc. The vertical dot-dashed line indicates the completeness limit at 2 kpc. Right bottom panel: Same as above 
for clumps with masses above the completeness limit at 2 kpc. 



the other properties, the mean values for the size, H2 vol- 
ume density, and surface density are 0.3 pc, 4.7x 10 5 cm" 3 , 
and 0.2 gem -2 for the clumps without MSX counterpart 
(the densities maybe also scaled by a factor ^2.5 like the 
mass, in case T d = 15 K), and 0.5 pc, 9.6 x 10 5 cm' 3 , and 
0.4 gem -2 for the clumps associated with mid-IR emis- 
sion. The median values are 0.3 pc, 4.0 x 10 4 cm" 3 , and 
0.11 gem -2 for the clumps without MSX counterpart, and 
0.4 pc, 3.6x 10 4 cm" 3 , and 0.14 gem" 2 for those with MSX 
counterpart. 

In order to assess whether the clumps without MSX 
emission are in an earlier evolutionary phase than the 
other ones, it would be necessary to map the cores in 
different molecular tracers to better derive their physi- 
cal properties, and in particular their temperature, which 
should be lower for precluster cores. 



5.4. The Mass Spectrum of the clumps 

Figure shows the histogram of the mass spectrum of 
the 1.2 mm clumps detected at a distance d < 6 kpc (top 
panels), and d < 2 kpc (bottom panels). The complete- 
ness limit, which has been estimated by calculating the 
mass of a 5cr detection at the upper distance limit and 
is shown as a vertical dot-dashed line in the left panels, 
is ~ 92 M© at 6 kpc, and ~ 10 M at 2 kpc. The right 
panels show the normalized cumulative mass distribution 
of clumps with masses above the completeness limit. The 
number of clumps above the completeness limit is 249 for 
d < 6 kpc, and 79 for d < 2 kpc. If the clump mass dis- 
tribution can be represented by a power law of the type 
dN/dM oc M~ a , then the histogram of the mass spec- 
trum can be fitted with a straight line of slope —a. The 
solid line in the figures corresponds to a — 2.35, i.e., the 
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Table 3. Power law indices of the mass spectrum of the clumps, dN/dM oc M~ a , for different mass ranges. The mass 
of the clumps has been derived either from the dust continuum or from gas emission. The arrows indicate the whole 
mass range for which the fit to the mass spectrum has been done. 
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Salpeter <|1955[l IMF; the dashed line to a = 1.7, cor- 
responding to the mass function of molecular clouds de- 
rived from gas, mainly CO, observations (e.g. Kramer et 
al. ©98); and the dotted line corresponds to the slope 
computed from a least squares fit to our data. This value 
is a = 2.10±0.02, with a fit correlation coefficient of 0.991, 
for the mass range 92 Af© <; Af <; 4100 Af© (d < 6 kpc), 
and a — 1.90 ± 0.04, with a fit correlation coefficient 
of 0.977, for the mass range 10 Af© <; M <; 387 Af© 
(d < 2 kpc). However, if ones does not take into account 
the last noisy point in the mass spectrum (M = 387 M©) 
when computing the least squares fit, and fits the clumps 
detected at d < 2 kpc with masses 10 Af© <; Af <; 120 Af©, 
the slope of the fitted mass spectrum is much flatter, 
a = 1.50, with a fit correlation coefficient of 0.999. Note 
that the errors reported for the slope values are those given 
by the fit, but do not represent the true uncertainties in 
the values of a, which depend also on other factors such 
as the binning used to determine the mass spectrum. 

As can be seen from Fig. which shows the mass 
spectrum and the best-fit power law for two different mass 
ranges, the mass spectrum, for masses above the complete- 
ness limit, cannot be fitted by a single power law for the 
whole mass range; that is, for 10 Af© <; Af < 4100 M©. 
The clump mass spectrum is slightly steeper for the high- 
mass end (M > 100 Af©), and it flattens for lower masses 
(M <; 100 Af©): for clumps detected at d < 6 kpc 
with masses above the completeness limit of 92 Af© (top 
left panel), the spectrum is well fitted with a — 2.1, 
consistent with a Salpeter-like power law; on the other 
hand, for clumps detected at d < 2 kpc, with masses 
10 Af© <; Af <; 120 M Q , the mass spectrum is better 
fitted with a = 1.5 or, if one takes also into account a 
clump with a mass M < 387 Af©, with a = 1.9, more 
consistent with the mass function of molecular gas clouds 
(e.g. Kramer et al. I1998f) . 

Such a change in the slope of the mass spectrum has 
been previously detected in other dust emission single-dish 
surveys of massive very YSOs (see Table |3J). Williams et 
al. I|2004fl have fitted the mass spectrum of their high- 
mass sample with a very flat slope, a — 1.14 for M < 
100 Mq, and with a = 2.32 for M > 100 Mq. Reid & 
Wilson (2005J have fitted the submillimeter clump mass 



function in NGC 7538 with a = 0.9±0.1 for the mass range 
15 Mq £ M £ 100 Mq, and with a = 2.0 ± 0.3 for the 
mass range 100 Af© <J M <^ 2700 Mq. In both cases, the 
breakpoint in the slope of the mass spectrum is well above 
the completeness limit of their samples, which suggests 
that the change in slope is real and not an observational 
effect. A power- law with a = 1.7 has been fitted by Tothill 
et al. H2002|) to the dust condensations with masses in the 
range 10 M© <J M <J 35 Mq in the Lagoon Nebula, which 
is the same region that we have mapped in our mosaic 
around 18014-2428 (Fig. |3J). Assuming that the break in 
the mass spectrum for massive dust clumps is real and not 
an artifact of the observations due to incompleteness, its 
shape is certainly interesting. 

For M J> 100 Mq, the slope of the mass spectrum 
of the dust clumps is a — 2.0-2.32 (see Table 0J). Such 
massive clumps are associated with pre- and protoclus- 
ters, so that the mass spectrum plotted in Fig. ^3 could 
correspond to the pre- and protocluster mass distribution. 
Interestingly, a slope of 2-2.32 is similar to the clump mass 
distribution in low-mass star-forming regions, which for 
masses 0.5 M© ^ M <; 10 M© mimics the stellar IMF 
(e.g. Testi & Sargent HUM Motte et al. ITM51 Johnstone 
et al. 120001 12001fl . The similarity between our result and 
those obtained in low-mass star-forming regions seems to 
suggest that the fragmentation of massive clumps may de- 
termine the IMF and the masses of the final stars. In other 
words, the processes that determine the clump mass spec- 
trum might be self-similar across a broad range of clump 
and parent cloud masses. Reid & Wilson (2005), based 
on theories of molecular cloud evolution (e.g. Gammie et 
al. EOOSl Tilley & Pudritz l2HM|) . suggest that turbulent 
fragmentation might be the dominant process that deter- 
mines the shape of the clump mass spectrum. 

For lower masses (10 Af© <; M <; 100 M©) the slope 
of the mass spectrum is shallower and consistent with the 
one found for clumps of similar masses observed in molec- 
ular lines (a ~ 1.7; Kramer et al. 11998(1 . This is also shal- 
lower than the spectrum of lower mass (Af <; 10 Af©) 
pre- and protostellar dust clumps. However, such a com- 
parison may be not appropriate because of the different 
tracer used in the line and continuum surveys, as line ob- 
servations are more sensitive to low-density material than 
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continuum imaging. Moreover, for our clumps the slope of 
the mass spectrum in the interval 10 Mq < M < 100 Mq 
is quite uncertain, because the mass range is relatively 
small and close to the completeness limit. Finally, the 
limited angular resolution of the single-dish observations 
of massive star-forming regions might be insufficient to 
resolve the clumps into the cores associated with indi- 
vidual (proto)stars. This makes it hazardous to compare 
the results with those of interferometric observations or 
with those of low-mass star-forming regions, which have 
a much higher angular resolution. In fact, Beuther & 
Schilkc (2004) resolved 12 clumps, with masses 1.7 Mq <; 
M <; 25 Mq, in their interferometric study of the massive 
star forming region IRAS 19410+2336 and determined a 
mass spectrum with a slope of 2.5, consistent with the 
stellar IMF. More studies of this type are needed to es- 
tablish the slope of the mass function in high-mass star 
forming regions. 

6. Conclusions 

We have extended to the southern hemisphere the project 
started by Palla et al. l|1991[l in the northern sky aimed 
at identifying high-mass protostellar candidates. We have 
carried out a 1.2 mm dust continuum emission survey with 
the bolometer array SIMBA at the SEST antenna of a 
sample of 235 sources selected from the IRAS-PSC: 93 
High sources with CS (Bronfman et al. I199fi[) . 125 Low 
with 6 < -30°, observed in CS and/or C 17 (with the 
exception of 15579-5347) by Fontani et al. ( 2005), plus 17 
Low with S > —30° observed in NH3 (with the exception 
of 18198-1429) by Molinari et al. ((TMB) . 

Massive dust clumps have been detected in all but 8 
regions, with usually more than one clump per region. 
The dust emission shows a variety of complex morpholo- 
gies, sometimes with multiple clumps forming filaments 
or clusters. Most of the sources with faint dust continuum 
emission are Low sources with no C 17 or CS detected 
towards the IRAS position. The mean clump has a lin- 
ear size of 0.5 pc, a mass of 320 Mq for a Td = 30 K, 
an H2 density of 9.5 x 10 5 cm -3 , and a surface density 
of 0.4 gcm~ 2 . The median values are 0.4 pc, 102 Mq, 
4 x 10 4 cm -3 , and 0.14 gem -2 , respectively. 

The mean value of the luminosity-to-mass ratio, L/M, 
is 99 Lq/Mq for all the sources, and 98 Lq/Mq and 
101 Lq/Mq for the Low and the High sub-samples, re- 
spectively. Such values are ~ 5 times lower than the aver- 
age value of known UC H11 regions with the same masses. 
This difference suggests that the sources in our sample, 
both Low and High, are in a younger pre-UC H11 phase, 
as the L/M ratio is expected to increase as the cores evolve 
and develop UC H11 regions. 

The physical properties of both Low and High sub- 
samples do not show significant differences (see Fig. 0). 
The luminosity of the High sources is on average higher 
than that of the Low sources (see Fig. |SJ), as are the lin- 
ear diameters and the masses of the clumps. The similar 



L/M ratio for both sub-samples suggests that they are in 
a similar evolutionary stage. 

There are 95 massive millimeter clumps in the sur- 
roundings of the candidate massive (proto)stars that are 
not associated with mid-IR MSX emission, either point- 
like or diffuse (see Fig- El - 70 out of the 95 clumps without 
MSX emission are clumps in maps around Low sources, 
which is the 16% of the Low clumps, while 20 are clumps in 
maps around High sources, 10% of the High clumps. Such 
clumps are potential prestellar or precluster cores, as one 
may expect that the bulk of their luminosity is emitted 
at millimeter and submillimeter wavelengths. The physi- 
cal properties of these clumps are similar to those of the 
MSX emitters, apart from the mass that is significantly 
lower than for clumps with MSX counterpart. However, 
such a difference could be due to the potential precluster 
clumps having a lower dust temperature. 

The mass spectrum of the clumps with masses above 
M ~ 100 Mq has been fitted with a power-law index 
a = 2.1, consistent with the Salpeter Ijl955|) stellar IMF, 
dN/dM ~ M~ 2 - 35 , or with the low -mass pre- and pro- 
tostellar dust clumps mass spectrum. On the other hand, 
the mass spectrum for clumps with masses 10 Mq <; M <J 
120 Mq is better fitted with a slope a — 1.5, more consis- 
tent with the mass function of molecular gas clouds (e.g. 
Kramer et al. H998|l . For M > 100 Mq, the massive dust 
clumps are probably tracing pre- or protoclusters, which 
have a mass spectrum similar to the stellar IMF. This sug- 
gests a self-similar process which determines the shape of 
the mass spectrum over a broad range of masses, from stel- 
lar to cluster size scales. For M <; 100 Mq, the shallower 
slope of the mass spectrum could be due to the limited an- 
gular resolution of single-dish observations, which is not 
enough to resolve the clumps into their real star-forming 
entities. 
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Table 1. Observed IRAS sources, alongside with source type, H(igh) or L(ow) according to the classification of Palla 
et al. 119911 position (precessed to J2000 coordinates), kinematic distance, and luminosity. The last column shows the 
sum of the mass estimated from the 1.2 mm continuum for clumps located at < 90" from the IRAS position. The 
luminosities and the masses have been estimated assuming the near distance, d ne ar, except where noted. 







a(J2000) 


<5(J2000) 


d a 

"'near 


^far 


L b 


M 


IRAS name 


Type 


h m s 


o / // 


(kpc) 


(kpc) 


(xlO 3 L ) 


(M e ) 


08140-3559 


H 


08 15 59.0 


-36 08 18 


3.8 


3.8 


10.6 


95 


08211-4158 


L 


08 22 52.3 


-42 07 57 


1.7 


1.7 


3.0 


74 
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L 


08 26 27.6 


-42 33 05 


1.4 


1.4 


1.5 


17 


08438-4340 


H 


08 45 36.0 


-43 51 01 
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4.3 
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08470-4243 
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15.2 
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08 49 26.7 


-43 17 13 
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c 
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3.2 
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(a) Kinematic distance estimated from the NH3 or CS line velocity using the rotation curve of Brand & Blitz (1993). 

(b) Luminosities calculated by integrating the IRAS flux densities. 

(c) No CS detected towards the IRAS source (Fontani et al. |5005). 

(d) No source detected at millimeter wavelengths. 

(e) Millimeter sources located at > 90" from the nominal IRAS point source position. 

(f) Far distance more than 200 pc from the galactic plane. 

(g) dnear < 100 pc. Luminosity and mass estimated assuming df ar . 

(h) No CS or C 17 observed. 

(i) Galactocentric distance < 2 kpc. 
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Table 2. Offset position, angular diameter, linear diameter, integrated flux density at 1.2 mm, mass, and density of 
the clumps detected towards each IRAS source. The offset positions are relative to the nominal IRAS point source 
positions. The last column indicates whether the clump is associated with MSX emission, either point-like or diffuse, or 
not. Linear diameters, masses, and densities of the clumps have been estimated assuming the near kinematic distance 
when the distance ambiguity could not be resolved, except where noted. 
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(a) The mass of the clump, M c i U mp, and the H2 volume density, nn 2 , are estimated from the 1.2 mm continuum for a dust 
temperature of 30 K, and using a gas-to-dust mass ratio of 100. 

(b) P: MSX point source within 40" from the millimeter clump peak emission position; D: diffuse MSX emission; N: No MSX 
emission. 

(c) No kinematic distance estimate. 

(d) No source detected at millimeter wavelengths. 

(e) MSX point source within 40" of nominal IRAS position, or diffuse MSX emission. 

(f) Unresolved source. 

(g) Linear size, mass and density estimated assuming the far kinematic distance. 

(h) The clumps listed have been detected in a field offset by 4.'7 from the nominal IRAS point source position. 
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Fig. 2. Overlay of the 1.2 mm continuum emission (contours) obtained with SIMBA at the SEST antenna, on the 
MSX emission at 21 fj,m (Band E) (image) towards the sources of our survey. The contour levels are 3<7, 6<r, and 
from 6cr to the maximum integrated emission in steps of 6er, where a is the rms noise of each map. The equatorial 
coordinates are in J2000 epoch. The noisy edges of the SIMBA map are visible. The black cross marks the nominal 
IRAS point source position. For IRAS 16428—4109, the imaged field is offset by A'.7 from the nominal IRAS point 
source position. For IRAS 18014-2428/18018-2426, the two IRAS sources correspond to Mol 36 (IRAS 18014-2428; 
black cross) and Mol 37 (IRAS 18018-2426: black trianale) from Molinari et al. dISSBl. 



